Zurek-Kibble Mechanism for the Spontaneous Vortex Formation in Nb — Al/Al ox /Nb 
Josephson Tunnel Junctions: New Theory and Experiment^ 
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New scaling behavior has been both predicted and observed in the spontaneous production of flux- 
ons in quenched Nb — Al/Alox/Nb annular Josephson tunnel junctions as a function of the quench 
time, tq. The probability /i to trap a single defect during the N-S phase transition clearly follows 
an allometric dependence on tq with a scaling exponent a = 0.5, as predicted from the Zurek- 
Kibble mechanism for realistic JTJs formed by strongly coupled superconductors. This definitive 
experiment replaces one reported by us earlier, in which an idealised model was used that pre- 
dicted o = 0.25, commensurate with the then much poorer data. Our experiment remains the only 
condensed matter experiment to date to have measured a scaling exponent with any reliability. 
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The Zurek-Kibble (ZK) scenario Q, HQ for continu- 
ous phase transitions proposes that transitions take effect 
as fast as possible i.e. the domain structure initially re- 
flects the causal horizons. This proposal is amenable to 
direct testing for transitions whose domain boundaries 
carry topological charge. This is the case for Josephson 
Tunnel Junctions (JTJs), where the topological charge 
is a fluxon i.e. a supercurrent vortex carrying a single 
quantum of magnetic flux <J> = h/(2e) in the plane of 
the oxide layer between the two superconductors that 
make up the JTJ. In the case considered in this paper 
of annular JTJs, obtained by the superposition of two 
superconducting rings, it corresponds to a magnetic flux 
that threads just one of the two rings an odd number of 
times. 

In 2000 an idealised model was proposed Q to test 
the ZK scenario for JTJs, assuming that causal horizons 
are constrained by the velocity of electromagnetic waves 
in the JTJ, the Swihart velocity 0, @- As a result, the 
probability fx for spontaneously producing one fluxon 
in the thermal quench of a symmetric annular JTJ of 
circumference C was predicted to scale with the quench 
time tq (the inverse quench rate) as 
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In Eq.QJ, £ is the Zurek-Kibble causal length, the cor- 
relation length of the relative phase angle at the time of 
defect formation. It is defined through Eq. (JTJ in terms of 
the cold correlation length £o, the relaxation time To of 
the long wavelength modes and the quench time tq, in 
turn defined by: T c /t q = -{dT /dt) T=Tc . Eq.Q holds 
for C < £. 



Under the assumptions of a) weak coupling of the super- 
conductors and b) exact critical slowing of the Swihart 
velocity at the critical temperature T — T c , we predicted 
a = 0.25 0. 

In 2001 our first proof-of-principle experiment was 
performed, to test the scaling law of Eq.Q. The exper- 
iment consisted of taking an annular JTJ isolated from 
its surroundings and making it undergo a forced phase 
transition by heating it above its superconducting criti- 
cal temperature and letting it to cool passively back to- 
wards the liquid He temperature. Once the thermal cy- 
cle is over, the junction I-V curve is inspected and any 
trapped fluxon can be detected by the appearance of cur- 
rent peaks at discrete voltages in the I-V characteristic of 
a JTJ. By detecting the voltage position of these peaks 
the number of trapped fluxons (or antifluxons) can be 
determined. 

The experiment of Ref. was remarkabl y su ccessful, as 
the onl y ex periment of the several |s 0, Hfl Ull lid 
IbH Hoi Till Iril previously performed to test scaling on 
condensed matter systems which was sensitive enough to 
show genuine ZK scaling behavior. However, this pre- 
liminary experiment suffered from severe restrictions; in 
particular a limited range of cooling rates, statistically 
poor data and above all insufficient shielding of the earth 
magnetic field with the possibility of systematic errors. 
The outcome was a commensurate with 0.25. 

In this paper we shall present results from a new ex- 
periment designed to circumvent these problems, that we 
shall discuss later. Our experiment shows extremely reli- 
able scaling behavior of the form Q), but with a = 0.5 to 
high accuracy. This is obviously at total variance with 
our earlier prediction. However, as we shall show, the 
value a — 0.25 depended critically on both assumptions 
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a) and b) cited earlier being exactly satisfied. For realistic 
JTJs these assumptions are only approximately satisfied 
in the immediate vicinity of T = T c . Since defects form 
very close to T c , we shall see that this is sufficient for the 
value of a to jump from 0.25 to a — 0.5 exactly. Thus, 
rather than its negation, our experiment provides an even 
more robust demonstration of the ZK scenario. 

To reiterate, the theory in Ref.Q was developed 
for JTJs whose electrodes are weak coupling supercon- 
ductors; for such JTJs the temperature dependence 
of the critical current density J A T) is given by the 
Ambegaokar-Baratoff equation : 



where 



7tA(T) A(T) 



(2) 



where A(T) is the superconducting gap energy and pn 
is JT J normal resistance area squared. Eq. J2J provides a 
linear decrease of J c near T c : 



■/.<T(/)l = n./,(l»(l-^ 



aJ c (0) — , (3) 



where T(t = 0) = T c and the dimensionless quantity 
a is approximately equal to 2A(0)/fc,eTc = 3.5. How- 
ever, our JTJs are based on Nb, a strong- coupling su- 
perconductor, for which Eq.(|2J) is not necessarily valid. 
In practice, high quality and reproducible barriers are 
achieved by depositing a thin Al overlay onto the Nb base 
electrode which will be only partially oxidized, leaving a 
Nb — Al bilayer underneath having a non-BCS tempera- 
ture dependence of the energy gap and of the density of 
states . The proximity effect in Nb - Al/Al ox /Nb JTJs 
has been extensively studied and it is known to influ- 
ence the electrical properties of the junctions, such as 
the current-voltage characteristic and the temperature 
dependence of the critical current density. Specifically, 
the proximity effect in S — N — I — S junctions is respon- 
sible for a subdominant temperature dependence of the 
critical current density 

HI 

in the vicinity of T c of the 

form: 



J c (T(t)) ~ a' J c (0) (l - Q =a'J c (0)(j-) , 
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where a' is a constant depending intricately on the degree 
of proximity. The last equation models the tail shaped 
dependence of J c vs. T near T c ; it has been theoretically 
derived and experimentally confirmed by Golubov et al. 
[20| in 1995. Rephrasing the arguments of Refill with 
Eq. |@} replacing Eq. © , the Josephson penetration depth 
Aj(T), which plays the role of system equilibrium coher- 
ence length £(T), near T c diverges linearly with time: 



Z{T{t))=\j{T{i)) 
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d s being electrode thickness . 

= d£(t)/dt < measures the rate at which the 
defects contract, i.e., the speed of interfaces between or- 
dered and disordered ground states. Since £(t) decreases 
with time t > 0, the earliest possible time t at which 
defects could possibly appear is determined by causality, 



m = -m, 
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where i is the causal time and c is the Swihart velocity. 
As we said, in Ref.0] we had assumed that the Swihart 
velocity vanished at T c whereas, for realistic JTJs, it just 
becomes very small. Swihart |5j has demonstrated that 
for a thin-film superconducting strip transmission line 
the solution for the velocity varies continuously as one 
passes through the critical temperature into the normal 
state. As a result, we assume c(t) = c nn near the tran- 
sition temperature where c nn is the speed of light in a 
microstrip line made of normal metals. In the case of a 
microstrip line made by two electrodes having the same 
thickness d s and the same skin depth 8, with d s << S, 
separated by a dielectric layer of thickness d ox and di- 
electric constant e, c nn ~ (2/ 5)^Jd ox d a / ep |5j. Its value 
depends on the temperature very weakly, but depends 
on the frequency / through 5 = yj p/irpf, p being the 
normal metal residual resistivity. 

The match c(T c ) = c nn is certainly realistic and we still 
have approximate critical slowing down insofar as c nn is 
much smaller than the zero temperature Swihart veloc- 
ity Co — \J d ox /2XLoep, i.e., when the zero temperature 
London penetration depth Alo << S 2 /d s . For 300nm 
thick Nb electrodes (p — 3.8pilcm and Alo = 90nm), 
d ~ 1mm at say / = 10kHz, so the last inequality is 
fully satisfied. At the same frequency, for a value of the 
specific barrier capacitance c s — e/d ox = 0.027-F/m 2 typ- 
ical of low current density Nb — Al/Al ox /Nb JTJs, we get 
Cnn = 6 • 10 3 m/s and c = 1.6 • 10 7 m/s |2l1 |. 
The solution of the causality equation Eq. (7J with a non- 
vanishing Swihart velocity yields: t = \f £,$tq / c n n- 
Inserting the value of t into Eq.|J5J we obtain the new 
Zurek length £: 



1/2 



(8) 



(5) 



where tq = £,o/c n n (to = 0(lns)). We reach the impor- 
tant conclusion for realistic JTJs that the probability /i 
for spontaneously producing one fluxon in the quench is 
still predicted to scale with the quench time tq accord- 
ing to Eq.JQ), but the critical exponent is now a = 0.5, 
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rather than a = 0.25. Detailed calculations will be re- 
ported elsewhere j^- 

However, it is worthy to observe that by varying tq in the 
experimentally achieved four decade range 1ms < tq < 
10s, we get 10^s < t < 1ms that is always much larger 
than To; it means that by the time the Josephson phase 
freezes the Josephson effect is well established. Further, 
in the same tq range the normalized freezing temperature 
T/T c at which the defect is formed is 0.99 < T/T c < 
0.9999. It would be really hard, if not impossible, to 
measure the temperature dependence of J c and c so close 
to Tq- We have then to resort to theoretical predictions. 

The new experiment has a faster and more reliable sin- 
gle heating system, obtained by integrating a meander 
line 50/im wide, 200nm thick, and 8.3mm long Mo resis- 
tive film in either ends of the 4.2mmx3mmx 0.35mm Si 
chip containing the Nb/AlOx/Nb trilayer JTJs. These 
resistive elements have a nominal resistance of 50f2 at 
LHe temperatures and, due to their good adhesion with 
the substrate, are very effective in dissipating heat in the 
chip. In fact, voltage pulses a few fis long and a few 
volts high applied across the integrated heater provided 
quench times as low as 1/zs, that is more than two or- 
ders of magnitude smaller than for the previous situation 
Q. Further, automatization of thermal cycles was imple- 
mented that allowed for much more robust statistics to be 
achieved. At the end of each thermal quench the junction 
I-V curve is automatically stored and an algorithm has 
been developed for the detection of the trapped fluxons. 
Finally, all the measurements have been carried out in a 
magnetic and electromagnetically shielded environment. 
The low frequency magnetic shielding was achieved by us- 
ing /z-metal, cryoperm and lead cans. During the quench 
the JTJ was also electrically isolated. 

The experimental results reported here are restricted 
to one of two annular JTJs with similar geometry to the 
sample used in the earlier experiment (C = 500/im and 
d s — 300nm) , but with about 50 times lower critical cur- 
rent density (J c (0) — 60A/cm 2 ): this leads to a value of 
£o = 17/xm (we have assumed a' = a = 3.5 in Eq.@). 
An extensive description of the chip layout, the experi- 
mental setup and the quenching data will be given else- 
where |22| while the details of the fabrication process can 
be found in Ref. |23j |. The quench time tq was contin- 
uously varied over more than four orders of magnitude 
(from 20s down to 1ms) by varying the helium exchange 
gas pressure inside the vacuum can and/or the width and 
the amplitude of the voltage pulse across the integrated 
resistive element. 

Fig. ^ shows on a log- log plot the measured frequency 
fi = n\/N of single fluxon trapping, obtained by quench- 
ing the sample N times for each value of a given quench- 
ing time tq, ri\ being the number of times that the 
inspection of the low temperature JTJ current-voltage 
characteristics at the very end of each thermal cycle 




FIG. 1: Log-log plot of the measured frequency /i of trap- 
ping single fluxons versus the quenching time tq . Each point 
corresponds to many thermal cycles. The vertical error bars 
gives the statistical error while the relative error bars in tq 
amounting to ±10% are as large as the dots' width. The 
solid line is the fit to an allometric relationship fi — aTg b 
which yields a = 0.0092 ± 10% (taking tq in seconds) and 
b = 0.51 ±6%. 



showed that one defect was spontaneously produced. N 
ranged between 100 and 2600 and ri\ was never smaller 
then 10, except for the rightmost point for which n\ = 3 
(and N = 1800). The sample has undergone a total 
of more than 100,000 thermal cycles without any mea- 
surable change of its electrical parameters. The vertical 
error bars gives the statistical error f\/y/n\. The mea- 
surement of tq by fitting solutions to the heat equation 
follows that of in its high accuracy. The relative error 
bars in tq amounting to ±10% are as large as the dot's 
width. 

To test Eq.©, we have fitted the data with an allomet- 
ric function fx = aTg 6 , with a and b as free fitting pa- 
rameters. A linear regression of log/i vs. log tq, rep- 
resented by the continuous line in Fig. ^ yields b — 
0.51 ± 6%, in excellent agreement with the predicted 
value 0.5. The same fit gives a = 0.0092 ± 10% (tak- 
ing tq in seconds) that is 6-7 times larger than the pre- 
dicted value C/VCoCnn (with C = 500/im, £o = 17//m 
and c nn = 6- 10 3 m/s). As a bound we only expect agree- 
ment in the overall normalization a to somewhat better 
than an order of magnitude. Empirically, the different 
condensed matter experiments have shown that the ratio 
^observed /flpredictcd varies widely from system to system; 
0(1) for superfluid 3 He H,0j very small for high-T c su- 
perconductors 0| . We point out that in our case the 
value of the prefactor a is dependent, although weakly, 
on the choice of a' and /, being C/ v / foBrm oc (a'//) 1 / 4 . 
The choice / = 10kHz was determined assuming that 
1/a 2 ~ 12kHz is the characteristic frequency of our sys- 
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tern at the time of the thermal quench. As we noted 
earlier, the new scaling exponent a = 0.5, obtained by 
applying causality arguments to realistic JTJs, is twice 
as large than that observed earlier 0, also in samples 
made with the same Nb — Al/Al ox /Nb technology. The 
reason for this discrepancy resides in the fact that at that 
time we were unaware of the high sensitivity of /i to ex- 
ternal magnetic fields which, although small, were most 
likely present. In the present experiment much care has 
been devoted both to avoid magnetic or current carry- 
ing materials in the cryoprobe during the quench and to 
screen the chip environment from 50Hz noise as well as 
from DC magnetic field. 

The data of Fig. ^ resolve another issue. It has been 
disputed that the Swihart velocity, which gives the be- 
havior above, is the relevant velocity for field ordering 
along the JTJ oxide. If we had taken the relevant velocity 
to be that of phase ordering in the individual supercon- 
ductors, as invoked by Zurek 0,0] when considering the 
spontaneous flux produced on quenching annuli of simple 
superconductors, at the same level of approximation we 
would have predicted b = 0.25 and the prefactor orders 
of magnitude larger. 

We have not observed the change in behavior predict- 
ing /i to change from the linear behavior with C/£ of 
Eq.Q to a random walk behavior in the phase, 



h 
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c 



-a/2 
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once /i is sufficiently large. A first guess would suggest 
that the transition from Eq.QJ to Eq.@ would occur 
when /i ~ 1. However, future experiments to be carried 
out on ring-shaped JTJS having the circumference larger 
than that of the sample reported in this paper and with 
tq in the same range of this experiment should clearly 
reveal the transition from Eq.Q to Eq.@. 

In summary, we see this experiment as providing un- 
ambiguous corroboration of Zurek-Kibble scaling over a 
wide range of quenching time tq in accord with our pre- 
dictions for Nb — Al/Al ox /Nb JTJs. As such, it replaces 
the experiment reported in Ref.Q by being more realis- 
tic theoretically and more sophisticated experimentally. 
We stress that this experiment is the only one to date 
to have confirmed the Zurek-Kibble scaling exponent for 
a condensed matter system. Further experiments can be 
devised to investigate the transition to the random walk 
regime and the effect of the thermal gradients. Such ex- 
periments are in the process of being performed. 
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